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ABSTRACT 

Kiloparsec-scale binary active galactic nuclei (AGNs) signal active supermassive black hole (SMBH) pairs 
in merging galaxies. Despite their significance, unambiguously confirmed cases remain scarce and most have 
been discovered serendipitously. In a previous systematic search, we optically identified four kpc-scale binary 
AGNs from candidates selected with double-peaked narrow emission lines at z = 0.1-0.2. Here we present 
Chandra and Hubble Space Telescope Wide Field Camera 3 (WFC3) imaging of these four systems. We 
critically examine and confirm the binary- AGN scenario for two of the four targets, by combining high angular 
resolution X-ray imaging spectroscopy with Chandra ACIS-S, better nuclear position constraints from WFC3 
F105W imaging, and direct starburst estimates from WFC3 F336W imaging; for the other two targets, the 
existing data are still consistent with the binary-AGN scenario, but we cannot rule out the possibility of only 
one AGN ionizing gas in both merging galaxies. We find tentative evidence for a systematically smaller X-ray- 
to-[Om] luminosity ratio and/or higher Compton-thick fraction in optically selected kpc-scale binary AGNs 
than in single AGNs, possibly caused by a higher nuclear gas column due to mergers and/or a viewing angle 
bias related to the double-peak narrow line selection. While our result lends some further support to the general 
approach of optically identifying kpc-scale binary AGNs, it also highlights the challenge and ambiguity of X- 
ray confirmation. 

Subject headings: black hole physics - galaxies: active - galaxies: interactions - galaxies: nuclei - galaxies: 
Seyfert - quasars: general - X-rays: galaxies 



1. INTRODUCTION 

1.1. Significance of Binary Supermassive Black Holes 

Most bulge-dominated galaxies harbor central supermas- 
sive black holes (SMBHs; Kormendy & Richstone 1995). As 
a result, binary 5 SMBHs are expected to form in galaxy merg- 
ers (Begelman et al. 1980; Milosavljevic & Merritt 2001; Yu 
2002); they are an inevitable and important consequence of 
the hierarchical bulge and SMBH formation process. Binary 
SMBHs are believed to have a significant dynamical impact 
on the nuclear stellar structure of massive elliptical galax- 
ies (e.g., Faber et al. 1997; Ravindranath et al. 2002; Gra- 
ham 2004; Merritt 2006; Kormendy & Bender 2009). The 
final inspiral and coalescence of hardened SMBH binaries 
are predicted to produce strong gravitational wave signals 
(Thorne & Braginskii 1976), the detection of which would 
offer a direct test of general relativity on cosmological scales 
(Thorne 1987). The identification and characterization of bi- 
nary SMBHs at various merger phases are valuable both for 
understanding galaxy/SMBH evolution and for probing fun- 



1 Based, in part, on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are associ- 
ated with program number GO 12363. 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, 
Cambridge, MA 02138 

3 Department of Astrophysical Sciences, Princeton University, Peyton 
Hall - Ivy Lane, Princeton, NJ 08544 

4 Einstein Fellow 

5 Following the initial nomenclature of Komossa et al. (2003) for NGC 
6240, we use "binary" AGNs to denote a pair of AGNs, also in line with the 
nomenclature "binary quasars" adopted in the literature. In this context, "bi- 
nary" does not necessarily presume that the black holes themselves are grav- 
itationally bound to each other (e.g., in the case of kpc-scale binary AGNs, 
the host galaxies dominate the potential well). 



damental physics (see a comprehensive review by Colpi & 
Dotti2011). 

1.2. Binary AGNs: the Kiloparsec Scales 

The frequency and statistical properties of binary active 
galactic nuclei (AGNs) may offer useful insights to the hier- 
archical merger paradigm of galaxy evolution (Yu et al. 201 1) 
and the role of mergers in AGN fueling 6 (Ellison et al. 201 1; 
Liu et al. 2012b). In a galaxy merger, if both BHs are simul- 
taneously accreting, they can be detected through spatially re- 
solved emission diagnostics which signal the presence of two 
AGNs. Theory suggests that merger-induced gas inflows be- 
come significant (therefore likely triggering AGN) at separa- 
tions under about a kpc (Hernquist 1989). Unlike the bound 
binary phase, which is still extremely challenging to image di- 
rectly (e.g., Burke-Spolaor 201 1), the "pairing" phase (where 
the separation between the two BHs, a, is a few tens pc to 
a few tens kpc) is the most accessible, because the two BHs 
are still resolvable at cosmological distances (typical separa- 
tion > 1"). Of particular interest is the late-pairing phase (a 
less than 10 kpc), which connects mergers in a cosmological 
context to pairs BHs in galaxies, and sets the stage for the 
subsequent evolution of close binaries. 



6 While the small-scale (~ 0.1-1 Mpc) quasar-quasar two-point correla- 
tion function suggests a clustering excess over the large-scale (> 1 Mpc) ex- 
trapolation (Hennawi et al. 2006; Myers et al. 2007; Hennawi et al. 2010; 
Shen et al. 2010), it is still unclear whether this is due to tidally enhanced 
BH accretion (e.g., Djorgovski 1991; Kochanek et al. 1999; Mortlock et al. 
1999), or is rather due to the small-scale clustering of their host dark matter 
halos (e.g., Hopkins et al. 2008; Richardson et al. 2012). The projected sep- 
arations of most of the observed binary quasars are on scales of tens of kpc 
and larger, which may still be too large for galaxy-galaxy tidal interactions 
to be effective (but see Green et al. 2010 for a counter-example of a 21 -kpc 
separation binary quasar observed to have tidal features indicative of ongoing 
interaction). 
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1.3. Systematic Searches for kpc-Scale Binary AGNs 

The past few years have seen significant increase in the in- 
ventory of kpc-scale binary AGNs, both from serendipitous 
discoveries and from systematic searches. While the exis- 
tence of kpc-scale binary AGNs has been confirmed in a few 
pioneering early discoveries (Komossa et al. 2003; Hudson 
et al. 2006; Bianchi et al. 2008) and further verified by more 
recent studies in X-rays (Brassington et al. 2007; Fabbiano 
et al. 201 1; Koss et al. 201 1; Mazzarella et al. 2012; see also 
Comerford et al. 201 1 for a candidate), radio (Fu et al. 201 lb), 
and optical broad emission lines (Shields et al. 2012), the fre- 
quency of occurrence and statistical properties of these sys- 
tems remain poorly constrained. 

Addressing the frequency and statistical properties of bi- 
nary AGNs requires systematic searches. A natural approach 
is to select candidates in galaxy mergers with double nuclei 
and follow up to identify binary AGNs using diagnostic obser- 
vations such as X-ray imaging (Guainazzi et al. 2005; Picon- 
celli et al. 2010; Koss et al. 2012; Teng et al. 2012) and/or spa- 
tially resolved optical spectroscopy (Barth et al. 2008; Com- 
erford et al. 2009b 7 ; Green et al. 2010; Greene et al. 2011; 
Liu et al. 201 1; Shields et al. 2012). However, because of the 
requirement that the two nuclei are resolved in ground-based 
optical imaging, systems identified using this approach are in 
general biased against the late-pairing phase. 

To mitigate this bias, an alternative approach is to se- 
lect candidates by kinematic signatures in spatially integrated 
spectra, in analogy to the case of spectroscopic binary stars. 
In particular, one such signature is the few hundred km s -1 
velocity splitting observed in AGN narrow emission lines 
(NELs; e.g., Sargent 1972; Heckman et al. 1981) such as 
[Om] AA4959,5007, which is seen in - 1% of low-redshift 
AGNs (Liu et al. 2010b; Smith et al. 2010; Wang et al. 2009; 
Ge et al. 2012). The working hypothesis is that the veloc- 
ity splitting signals the projected relative orbital motion of 
two narrow line regions (NLRs), each ionized by its own cen- 
tral AGN (e.g., Zhou et al. 2004; Gerke et al. 2007; Comer- 
ford et al. 2009a; Xu & Komossa 2009; Barrows et al. 2012). 
By selection, only binaries with projected angular separations 
smaller than the spectroscopic aperture size will be included. 
In principle, binary AGNs with separations as small as a few 
tens of pc (limited by the intrinsic size of NLRs, which are 
~ 50-1000 pc, and scale approximately as L 5 ; Bennert et al. 
2002) may be identified, if the associated double stellar nuclei 
are resolvable by followup observations using higher resolu- 
tion near infrared (NIR) imaging with HST and/or ground- 
based adaptive optics (AO). Therefore, this kinematics ap- 
proach should be well suited for identifying binary SMBHs 
in the late-pairing phase. 

A major obstacle in identifying binary AGNs using the nar- 
row emission-line splitting signature, however, is that such 
profiles can also arise from NLR gas kinematics around sin- 
gle AGNs, such as rotating disks or bi-conical outflows (Axon 
et al. 1998; Veilleux et al. 2001; Crenshaw et al. 2010; Shen 
et al. 2011; Rosario et al. 2010; Fischer et al. 2011; Smith 
et al. 2011, 2012). Nevertheless, our followup observations 
(Liu et al. 2010a; Shen et al. 2011) of a subset of a system- 
atically selected sample of 167 AGNs with double-peaked 
NELs (Liu et al. 2010b; see also Smith et al. 2010; Wang 
et al. 2009), as well as studies by other groups (McGurk et al. 

7 But see Civano et al. (2010, 2012) and Blecha et al. (2012a) for an alter- 
native explanation for this particular candidate. 



2011; Fu et al. 2012), have demonstrated the feasibility and 
importance of combining higher resolution NIR imaging (Fu 
et al. 2011a; Rosario et al. 2011) and spatially resolved opti- 
cal spectroscopy (Comerford et al. 2012) to discriminate kpc- 
scale binary AGNs from single-AGN-NLR gas kinematics. 
Kpc-scale binary AGNs tend to show two concentrated [O in] 
nuclei spatially coincident with two stellar bulges in a merger, 
with the dynamics dominated by the potential of the individ- 
ual stellar bulges. Objects with complex NLR kinematics in 
contrast, usually exhibit bi-cone/disk shaped diffuse [Om] 
with a smooth single-peaked stellar background, as well illus- 
trated by the local example Mrk 78 (Whittle & Wilson 2004; 
Fischer et al. 201 1). Roughly 10% of the objects we have fol- 
lowed up are best explained by binary AGNs at (projected) 
kpc-scale separations. 

1.4. This Work: X-ray Imaging of Optically Identified 
kpc-Scale Binary AGNs 

Here and in a companion paper (Liu et al. 2012a, hereafter 
Paper II), we present F336W/Z/- and F105W/Y-band images 
obtained using Wide Field Camera 3 (WFC3) on board the 
Hubble Space Telescope (HST), and 0.5-10 keV X-ray im- 
ages taken with the Chandra X-ray Observatory (Weisskopf 
et al. 1996) Advanced CCD Imaging Spectrometer (ACIS; 
Garmire et al. 2003) of the four optically selected kpc-scale 
binary AGNs identified by Liu et al. (2010a). Although 
our ground-based NIR imaging and spatially resolved optical 
spectroscopy strongly suggest that these galaxy mergers host 
binary AGNs (Liu et al. 2010a; Shen et al. 201 1), the case is 
not watertight. While seven of the eight narrow emission-line 
nuclei in the four galaxies are optically classified as Type 2 
Seyferts, one [Om] -faint nucleus is optically classified as ei- 
ther a Type 2 Seyfert, a LINER, or a LINER-H II composite, 
the latter two cases of which may be due to starburst and/or 
shock heating rather than AGN excitation (e.g., Lutz et al. 
1999; Terashima et al. 2000; Eracleous et al. 2002). More im- 
portantly, even in the "Seyfert-Seyfert" cases, there could be 
only one AGN, which ionizes gas both merging components. 
Photo-ionization arguments based on spatially resolved opti- 
cal spectroscopy were unable to rule out this possibility (Liu 
et al. 2010a), given that the nuclear separation is not much 
larger than the sizes of individual NLRs, and the large system- 
atic uncertainties in the electron density measurements from 
diagnostic emission line ratios. 

Our present work is motivated to further clarify these ambi- 
guities concerning the nature of the ionizing sources in opti- 
cally selected kpc-scale binary AGNs. HSTAVFC3's F-band 
images allow us to get better positional priors to resolve the 
closely separated double nuclei in X-ray imaging; the U- 
band images offer constraints on spatially resolved star for- 
mation activity in the host galaxy (Section 3.1). Utilizing 
Chandra ACIS's superb spatial imaging spectroscopy capa- 
bility in the X-rays (Section 3.2), we put more direct con- 
straints on the intrinsic X-ray luminosity for each individual 
nucleus in the merging galaxies than those estimates inferred 
empirically from [O in] A5007 emission-line luminosity (Sec- 
tion 4.1). Combined with constraints on the contribution from 
star formation to the observed X-ray luminosities estimated 
from U -band imaging (Section 4.2), we critically examine 
the purported binary- AGN nature of our targets (Section 4.3). 
Combining the new X-ray observations with our previous op- 
tical spectroscopy, we tentatively characterize the fraction of 
optical binary AGNs that are weak or Compton-thick X-ray 
emitters, and compare with the fraction among single AGNs 
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TABLE 1 

KPC-SCALE BINARY- AGN CANDIDATES IMAGED WITH HST/WFC3 AND Chandra X-RAY OBSERVATORY/ ACIS. 



Target Name 
(1) 


Plate 

(2) 


Fiber 

(3) 


MJD 

(4) 


Redshift 

Zc 
(5) 


HST 
Obs. UT 
(6) 


A6» Y 
(") 
(7) 


r,,.Y 
(kpc) 
(8) 


Chandra 
Obs. UT 
(9) 


A6> X -ray 

(") 
(10) 


SDSS J110851.04+065901.4 


1004 


182 


52723 


0.1816 


20110513 


0.70 


2.1 


20110210 


0.82 


SDSS J113126.08-020459.2 


327 


394 


52294 


0.1463 


20110524 


0.70 


1.8 


20110211 




SDSS J114642.47+511029.6 


881 


241 


52368 


0.1300 


20110620 


2.71 


6.3 


201 10423 


3.03 


SDSS J133226.34+060627.4 


1801 


250 


54156 


0.2070 


20110310 


1.50 


5.1 


20111128 





NOTE. — Col. (1): SDSS names with J2000 coordinates given in the form of "hhmmss.ss+ddmmss.s". Cols. (2)-(4): SDSS spectroscopic plate number, fiber ID, and Modified 
Julian Date. Col. (5): systemic redshift measured from stellar continuum absorption features in the SDSS fiber spectra. Cols. (6) & (9): dates of the HST and Chandra observations. 
Cols. (7) & (8): projected angular and physical separation between the double nuclei measured from HST 7 -band images. Col. (10): angular separation between the double nuclei 
measured from ACIS X-ray images. See Table 2 for measured positions of the individual nuclei. 



(Section 4.4). Since our targets were selected in a systematic 
search, our results have general implications for the general 
approach of identifying kpc-scale binaries in double-peaked 
AGNs (Section 5.1), the attributes and limitations of optical 
identification compared to X-ray searches (Section 5.2), and 
the frequency of kpc-scale binary AGNs (5.3). 

Throughout this paper, we assume a concordance cosmol- 
ogy with £l m = 0.3, Qa = 0.7, and Ho = 70 km s" 1 Mpc" 1 , and 
use the AB magnitude system (Oke 1974). 

2. TARGET SELECTION 

In Table 1 we list basic photometric and spectroscopic prop- 
erties of the four binary-AGN candidates. These candidates 
were discovered from a subset of a parent sample of 167 
Type 2 AGNs with double-peaked [Om] AA4959,5007 emis- 
sion lines (Liu et al. 2010b). The parent sample was iden- 
tified from 14,756 Type 2 AGNs optically selected from the 
spectroscopic catalog of the Sloan Digital Sky Survey (SDSS; 
York et al. 2000) Data Release Seven (DR7; Abazajian et al. 
2009). We focused on Type 2 (i.e., obscured) AGNs which 
allow us to study the host galaxy properties without much 
contamination from the AGNs. The optical emission line ra- 
tios [Om] A5007/H/3 and [Nil] A6584/Ha are characteristic 
of Type 2 Seyferts according to the Kewley et al. (2001) cri- 
terion based on the Baldwin et al. (1981) diagnostics. We 
conducted ground-based deep near-infrared (NIR) images and 
optical slit spectra from the Magellan 6.5 m and the Apache 
Point Observatory 3.5 m telescopes (Liu et al. 2010a; Shen 
et al. 201 1). We identified four strong kpc-scale binary AGN 
candidates out of 43 objects observed 8 (Liu et al. 2010a). In 
each system, the NIR images reveal tidal features and dou- 
ble stellar components with a projected separation of several 
kpc, while optical slit spectra show two Type 2 Seyfert nuclei 
(except for one nucleus which could also be a LINER or a 
composite) spatially coincident with the stellar components, 
with line-of-sight velocity offsets of a few hundred km s" 1 . In 
Table 2 we list redshift and [Om] A5007 (hereafter [Om]) lu- 
minosity measurements for each individual nucleus from our 
ground-based longslit spectroscopy. 

3. OBSERVATIONS, DATA REDUCTION, AND ANALYSIS 

3.1. HST/WFC3 F336W and F105W Imaging 

The four optically selected kpc-scale binary AGNs were ob- 
served using the WFC3 on board the HST in Cycle 18 (Pro- 
gram: GO 12363, PI: Shen). Each target was imaged in the 
UVIS/F336W (U band; Dressel 2010) and IR/F105W (wide F 

8 See also a fifth candidate, SDSS J1356+1026, reported by Shen et al. 
(201 1) and Greene et al. (201 1, 2012). 



band) filters within a single HST orbit. We refer to Paper II for 
details of our HST observations. Here we briefly describe the 
data relevant for addressing the nature of the ionizing sources. 

The U- and F-band images were calibrated both photomet- 
rically and astrometrically. The typical relative astrometric 
accuracy is 0."004 for the f/-band and 0."01 for the F-band 
images. To improve absolute astrometric accuracy and to 
compare with X-ray images, we have registered the U- and 
F-band images with the SDSS astrometry. The resulting abso- 
lute astrometric uncertainties of the registered U - and F-band 
images were estimated as ~ 0."25 and <~ 0."20, respectively 
(Paper II). We list the F-band nuclear positions in Table 2 and 
the inferred separations between the double nuclei for each 
target in Table 1 . These F-band nuclear separations agree with 
those measured from our ground-based NIR imaging within 
uncertainties. 

At the redshifts of our targets (0.130-0.207; Table 1), the 
U -band filter covers rest-frame <~ 2600-3200 A. For obscured 
AGNs, this wavelength range is dominated by continuum 
photospheric emission from host galaxy young stellar popu- 
lations, which offers a useful indicator for star formation rate 
(SFR; e.g., Cram et al. 1998). The U -band filter also covers 
line emission from ionized gas, although the contamination 
from even the strongest line, Mg n A2800, is likely insignif- 
icant (<1%), given that its typical equivalent width is small 
(e.g., 5.2 ± 0.8 A, as measured from the composite spectrum 
of Type 2 AGNs by Zakamska et al. 2003). The U -band im- 
ages may also contain AGN light from the obscured nuclei 
scattered into our line-of-sight by dust and/or gas (e.g., Za- 
kamska et al. 2006), but our data suggest that the contribu- 
tion is likely insignificant (<5%), because of the moderate 
AGN luminosities of our sample (Table 2) and the absence of 
a broad H/3 or Ha component in the optical spectra (e.g., Liu 
etal. 2009). 

Our high resolution U -band imaging provides constraints 
on the intensity and spatial distribution of star formation ac- 
tivity in the host galaxies (Paper II). The inferred nuclear 
SFRs are useful for estimating the X-ray contribution from 
star formation-related processes (Section 4.2). In table 3 we 
list the observed f/-band flux fu and luminosity Lff s inte- 
grated within a region centered on each F-band nucleus posi- 
tion, with an aperture size matched to that of X-ray extraction 
(Sections 3.2.2 and 3.2.3). To estimate the intrinsic L v , we 
have adopted an extinction correction for each nucleus (Ta- 
ble 3), based on the Balmer decrement F-n a /F-np measured 
from our spatially resolved optical spectroscopy (Liu et al. 
2010a). We assume the extinction curve of Cardelli et al. 
(1989) with R v = 3.1 to calculate the Balmer decrement. We 
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FIG. 1 . — Chandra ACIS X-ray images (0."25 pixel -1 binning; unsmoothed) of the four optically selected kpc-scale binary AGNs. The left, middle, and right 
columns show the full band, zoomed-in soft and hard bands, respectively. Contours are HST/WFC3 F105W (Y-band; in red) and F336W ((/-band; in green) 
images in linear spacing. Black crosses indicate F-band nuclear positions. Typical absolute (relative) astrometric uncertainty is 0."20 (0."15) for our ACIS 
images (see Section 3.2.1 for details) and 0."20 (0."01) for the 7-band images. Major tickmarks are separated by l."0. 
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have carefully subtracted the host galaxy stellar continuum 
using population synthesis models (Bruzual & Chariot 2003) 
with the fitting method of Liu et al. (2009), to avoid the ratio 
Fyia/Fyip being overestimated due to strong Balmer absorp- 
tion in post-starburst populations. To estimate SFR from L v , 
we have adopted the empirical calibration of Hopkins et al. 
(2003), which is given by 



Moyr" 1 V 1-81 x 10 28 ergs s-'Hz" 1 / 

with an rms scatter of 0.13 dex. The relation was based on 
the SDSS M-band luminosity (corrected for obscuration) of 
2625 star-forming galaxies in the SDSS DR1 (Abazajian et al. 
2003), which was calibrated against SFR estimates inferred 
from Ha emission-line measurements (corrected for aperture 
and obscuration effects) according to the relation of Kenni- 
cutt (1998); it is valid 9 for 2 x 10 28 ergs s^Hz -1 < L v < 
10 30 ergs s _1 Hz _1 . Hopkins et al. (2003) have shown that the 
inferred SFR V is consistent with SFR estimates from the 1 .4 
GHz luminosity (which was in turn calibrated from the FIR 
luminosity according to Kennicutt 1998; Bell 2003) with an 
rms scatter of 0.23 dex. Below in Section 4.2.1 we discuss 
systematics and uncertainties in our SFR estimates and by ex- 
tension the inferred X-ray luminosity due to star formation- 
related processes (Section 4.2). 

3.2. Chandra ACIS X-ray Imaging Spectroscopy 

The four kpc-scale binary AGN candidates were observed 
with the ACIS-S on board the Chandra X-ray Observatory 
between 2011 February and November (Cycle 12 Program: 
G01-12127X, PI: Shen). Exposure times ranged from 16 
ks to 24 ks (Table 4). They were set by the requirement of 
obtaining ~ 100 counts in the 0.5-10 keV from the weaker 
nucleus of each target. The counts were estimated from the 
[O in] luminosity for each nucleus, using the empirical corre- 
lation between 2-10 keV (unabsorbed) and [Om] luminosi- 
ties (extinction corrected) from Panessa et al. (2006) as the 
baseline value, taking into account systematic uncertainties 
using the Heckman et al. (2005) relation (for observed lumi- 
nosities) for optically selected single Type 2 AGNs, assum- 
ing a single power-law spectrum with an absorbing column 
density Nh = 10 23 cm" 2 (typical for Type 2 Seyferts; Bassani 
et al. 1999) and a photon index T = 1 .8 (typical for unabsorbed 
Seyferts Green et al. 2009). All the targets were observed on- 
axis on the S3 chip, 9" to 15" away from the aimpoint. We 
examined the light curves and found no flares of either the 
sources or the background in each of the observations. 

We reprocessed the data using the standard Chandra Inter- 
active Analysis of Observations (CIAO) software (Fruscione 
et al. 2006) with version 4.4. We ran the chandra_repro 
script on the standard Level 2 event file for the recommended 
processing steps by the Chandra X-ray Center, applying the 
latest calibration (CALDB 4.4.1). The process corrected for 
charge transfer inefficiency and time-dependent gain. The 
energy-dependent subpixel event repositioning (EDSER) al- 
gorithm (Li et al. 2004) was applied to improve the image 
quality of ACIS-S data for sources near the optical axis of the 

9 For the three weakest {/-band nuclei in our sample (Table 3), the Lu 
measurements are ~ 10 times lower than the faint luminosity end of the SDSS 
star-forming galaxy sample studied by Hopkins et al. (2003). To estimate 
SFRs for these faint nuclei, we assume the SFR-L[/ relation by Hopkins et al. 
(2003) extrapolated to lower luminosities. 



telescope, where the PSF is under-sampled by the 0."5 ACIS 
pixels. Our targets have low count rates, so pileup effects 
(Ballet 1999) were insignificant. 

3.2.1. Astrometric Uncertainty of ACIS Images 

First we discuss astrometry of ACIS images and our effort 
of obtaining accurate alignment between the X-ray and opti- 
cal images, which are important to determining the nature of 
X-ray sources. For sources within 3 arcmin of the aimpoint, 
the typical absolute ACIS-S astrometric accuracy is < 0."6 10 
(radius size of the overall 90% uncertainty circle of ACIS-S 
absolute position; Chandra Proposers' Observatory Guide 11 , 
hereafter POG). The relative astrometric accuracy is 0."15 
(90% limit) for on-axis sources 12 . 

Astrometric calibration was applied as part of the pipeline 
processing of ACIS images. To verify the astrometric accu- 
racy, we ran wavdetect (Freeman et al. 2002), which is a 
wavelet-based algorithm for spatial analysis of Poisson data, 
to detect sources as references for any fine alignment, if 
needed, between the ACIS images and the SDSS, to which 
our F-band images have been registered. We used mkpsfmap 
to create observation-specific PSF map files instead of us- 
ing the PSF table. We adopted a high significance threshold 
(sigthresh=10~ 8 , corresponding to one spurious source in a 
10 4 x 10 4 pixel map) to ensure robust source detection. 

For SDSS Jl 108+0659 (SDSS Jl 131-0204), the astrome- 
try of the ACIS image agrees with that of the SDSS within 
0."2 (0."2), based on three (six) SDSS-matched sources de- 
tected by Chandra within 2 arcmin (5 arcmin) of the aimpoint. 
For SDSS Jl 146+51 10, only one bright X-ray source was de- 
tected in the FOV, 1 .6 arcmin away from the aimpoint, whose 
ACIS and SDSS positions agree within 0."2. For SDSS 
J1332+0606, the astrometry of the ACIS image agrees with 
that of the SDSS within 0."3, based on four SDSS-matched 
sources detected within 3 arcmin of the aimpoint. Given these 
results, we do not apply any further astrometry correction for 
the ACIS images, because the agreement is already compara- 
ble to the SDSS astrometric accuracy. 

Figure 1 shows the unsmoothed ACIS images of our four 
targets in the full (0.5-10 keV), soft (0.5-2 keV), and hard 
(2-10 keV) bands, respectively. Given the low count levels 
(Table 4), we do not apply any smoothing to avoid artifacts. 
As we will show in Section 4.4, our targets are significantly 
weaker hard X-ray emitters than those predicted from both 
the Panessa et al. (2006) and Heckman et al. (2005) relations 
for single optically selected AGNs (by ~ 0.8 ± 0.2 dex and 
~ 1 .9 ± 0.3-2.4 ± 0.3 dex in observed and unabsorbed 2-10 
keV luminosities, respectively), resulting in far fewer counts 
than we expected. Five of the eight nuclei in our targets were 
detected in the full band (both nuclei in SDSS Jl 108+0659, 
both nuclei in SDSS Jl 146+51 10, and the NE nucleus in 
SDSS J 1332+0606), of which four were detected in both soft 
and hard bands (both nuclei in SDSS Jl 108+0659, the SW 
nucleus in SDSS Jl 146+51 10, and the NE nucleus in SDSS 
J 1332+0606), whereas one was only detected in the soft band 
(the NE nucleus in SDSS J 1 146+5 110). The other three nuclei 

10 This was inferred based on measuring the distances between the Chan- 
dra X-ray source positions and corresponding otical/radio counterpart posi- 
tions from the Tycho2 (with astrometric accuracy of ~ 25 mas; H0g et al. 
2000) and ICRS (with astrometric accuracy of ~ 1 mas; Ma et al. 1998) cat- 
alogs. 

11 http://cxc.harvard.edu/proposer/POG/; see also 
http://cxc.harvard.edu/cal/ASPECT/celmon/. 

12 based on the 900 ksec ACIS-I observation of the Orion Nebula; POG. 
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TABLE 2 

HST/WFC3 F-band and Chandra/ACIS X-ray 0.2-10 keV Astrometry of the Double Nuclei 





Kedsniit 


i _ _ j obs 
lo g L [OIII] 


lo g L [OIII] 


R.A.y 


Dec.y 


R-A.X-ray 


DeC.X-ray 




Object Name 


Ze 


(ergs s ') 


(ergs s ') 


(J2000) 


(J2000) 


(J2000) 


(J2000) 


(") 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SDSS J1108+0659NW 


0.1812 


42.16 


42.66 


11:08:51.029 


+06:59:01.32 


11:08:51.031 


+06:59:01.26 


0.07 


SDSS J1108+0659SE 


0.1820 


41.52 


41.64 


11:08:51.061 


+06:59:00.81 


11:08:51.069 


+06:59:00.66 


0.19 


SDSS J1131-0204W 


0.1454 


41.40 


42.34 


11:31:26.042 


-02:04:59.33 








SDSS J1131-0204E 


0.1470 


41.31 


42.22 


11:31:26.088 


-02:04:59.21 








SDSS J1146+5110SW 


0.1293 


41.93 


42.18 


11:46:42.466 


+51:10:29.46 


11:46:42.504 


+51:10:29.45 


0.36 


SDSS Jl 146+51 10NE 


0.1303 


41.38 


41.65 


11:46:42.630 


+51:10:31.69 


11:46:42.672 


+51:10:32.03 


0.52 


SDSS J1332+0606SW 


0.2057 


41.06 


41.64 


13:32:26.340 


+06:06:27.31 








SDSS J1332+0606NE 


0.2074 


41.83 


42.84 


13:32:26.372 


+06:06:28.73 


13:32:26.364 


+06:06:28.63 


0.16 



NOTE. — Col. (2): emission-line redshift measured from spatially resolved optical spectra (Liu et al. 2010a). Col. (3): observed [O III] A5007 emission-line luminosity measured 
from spatially resolved optical spectra (Liu et al. 2010a). Cols. (4) & (5): coordinates of the double nuclei measured from HST F-band images. Typical absolute (relative) astrometric 
uncertainty is 0."2 (0."01); Cols. (6) & (7): coordinates measured from ACIS images. Typical absolute (relative) astrometric uncertainty is 0."2 (0." 15). See Section 3.2.1 for details; 
Cols. (8) & (9): difference between the X-ray and F-band measured positions. 



TABLE 3 

Nuclear {/-band Fluxes and Luminosities, Extinction Estimates, and Inferred Star Formation Properties 





log/S bs 


logL° bs 


Fh/3 

(4) 


E(B-V) 


Au 


logLfv 


SFRu 


L SF 

^X,0.5-2 keV 


L SF 
(10 40 ergs 


Object Name 


(Jy) 


(ergs s -1 Hz -1 ) 


(mag) 


(mag) 


(ergs s" 1 Hz -1 ) 


(M yr 1 ) 


(10 40 ergss-') 


(1) 


(2) 


(3) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


SDSS J1108+0659NW 


-4.55 


28.4 


4.1 


0.4 


1.9±0.7 


29.2±0.3 


10 


6 


6 


SDSS J1108+0659SE 


<-5.93 


<27.0 


3.1 


0.1 


0.5±0.8 


<27.2±0.3 


<0.1 


<0.03 


<0.03 


SDSS J1131-0204W 


<-5.97 


<26.8 


5.6 


0.7 


3.8±0.3 


<28.3±0.1 


<1 


<0.5 


<0.6 


SDSS J1131-0204E 


-5.69 


27.1 


5.5 


0.7 


3.7±0.4 


28.5±0.2 


2 


1 


1 


SDSS J1146+5110SW 


-4.91 


27.7 


3.1 


0.1 


0.4±0.1 


27.9±0.1 


0.4 


0.2 


0.2 


SDSS Jl 146+51 10NE 


-4.44 


28.2 


3.1 


0.1 


0.5±0.1 


28.4±0.1 


2 


0.7 


0.8 


SDSS J1332+0606SW 


-5.39 


27.7 


3.5 


0.2 


0.9±0.3 


28.1 ±0.1 


0.6 


0.3 


0.3 


SDSS J1332+0606NE 


-5.25 


27.8 


4.7 


0.5 


2.3±0.9 


28.8±0.4 


4 


2 


2 



NOTE. — All measurements for each nucleus were made within the same aperture as for the X-ray extraction. Col. (2): observed U -band flux density. Col. (3): observed (7 -band 
luminosity density. Col. (4): Balmer decrement measured from our ground-based slit spectra (Liu et al. 2010a). Col. (5): color excess estimated from the Balmer decrement. Col. 
(6): f/-band extinction estimated from the Balmer decrement and uncertainty due to aperture coverage mismatch (see Section 4.2.1). Col. (7): extinction-corrected U -band luminosity 
density and uncertainty (propagated from the extinction uncertainty due to aperture mismatch only, not including that due to uncertain dust geometry). Col. (8): star formation rate 
inferred from the extinction-corrected (7-band luminosity density (see Section 3.1 for details). Typical uncertainty is ~ 0.3 dex, which was estimated by convolving that propagated 
from Lu with the rms scatter of the SFR-Lu calibration (Equation 1). Cols. (9) & (10): X-ray luminosities inferred from the SFR estimate, assuming the empirical calibration of Ranalli 
et al. (2003, see Section 4.2 for details). Typical uncertainty is ~ 0.4 dex, which was estimated by convolving that propagated from SFR estimates with the rms scatter of the SFR-Lx 
calibrations (Equations 3 and 4). 



were undetected in the X-rays (both nuclei in SDSS J 1 131- 
0204, and the SW nucleus in SDSS J1332+0606). 

Of the three X-ray detected targets, the F-band nuclei 
in SDSS Jl 146+51 10 and in SDSS J1332+0606 are sepa- 
rated by l."5 and 2. "7, respectively (Table 1), which are 
well within the resolving power of ACIS; the F-band nu- 
clei in SDSS Jl 108+0659 require more careful decomposi- 
tion (Section 3.2.3). For the X-ray detected nuclei in SDSS 
Jl 146+51 10 and SDSS J1332+0606, we compare their po- 
sitions with the F-band nuclear positions. We measure the 
X-ray position with a 2D image fitting analysis using Sherpa 
(Freeman et al. 2001). We adopt a constant for the background 
and the PSF images as convolution kernels to fit the sources. 
We applied the F-band positions as the initial guesses. PSF 
images were created with the MARX software 13 using the 
PSF-ray table generated by the Chandra Ray Tracer (ChaRT; 
Carter et al. 2003), which simulates the best available PSF at 
any off-axis angle and for any energy or spectrum. As listed 
in Table 2, the X-ray positions agree with the F-band posi- 
tions within the uncertainties for all the detected sources. We 
have compared the radial profile of our targets against PSF 

13 http://space.mit.edu/ASC/MARX/index.html 



models. Each nucleus component is consistent with being an 
unresolved point source. 

3.2.2. Source Extraction for Well Separated Nuclei 

For SDSS Jl 146+51 10 and SDSS J1332+0606, we run 
dmextract to extract X-ray counts for each individual nucleus. 
We use the positions and their uncertainties measured for each 
nucleus from the HST F-band images as priors for source ex- 
traction. The background counts were extracted from source- 
free annular regions concentric to the target regions. We re- 
port the number of background subtracted counts in separate 
soft (0.5-2 keV) and hard (2-10 keV) band and in the full 
band in Table 4. 

For each of the two nuclei in SDSS Jl 146+51 10, we 
adopted a circular region with a 1 ."2 radius for source extrac- 
tion. The adopted size ensured no overlap between the two 
sources, with each region containing more than 95% of the 
encircled energy (two dimensional integral of the PSF). For 
SDSS J 1332+0606, we extracted the counts for each nucleus 
in circular regions of l."0 radii, to ensure that there was no 
overlap between the extraction regions of the two nuclei. No 
significant detection was obtained for the southern nucleus in 
SDSS J 1332+0606; we estimated a 3-er upper limit accord- 



Kpc-Scale Binary AGN I. Nature of the Ionizing Sources 



7 



TABLE 4 

X-ray Properties and Spectral Models 





N a (Galactic) 


Exposure 


T 


5 


H 






A/h(T = 1.8) 


Object Name 


(10 20 cm" 2 ) 


(seconds) 


(counts) 


(counts) 


(counts) 


HR 


r 


(10 22 cm" 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 



SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 



J1108+0659NW 

J1108+0659SE 

J1131-0204W 

J1131-0204E 

J1146+5110SW 

J1146+5110NE 

J1332+0606SW 

J1332+0606NE 



3.99 
3.99 
3.34 
3.34 
1.58 
1.58 
2.21 
2.21 



19273 
19273 
23723 
23723 
15432 
15432 
23952 
23952 



25.4±5.0 
12.4±3.5 
<6.6 
<6.6 
16.6±5.2 
4.6±3.4 
<10.9 
7.7±2.9 



20.5±4.5 
2.5±1.6 



2.9±2.9 
4.6±3.4 

3.9±3.2 



4.3±2.1 
9.1±3.0 



13.6±4.8 
<1.8 

3.7±3.2 



-0.60^" 
0.46^ 



0.65*;!? 
<-0.81 

-0 (TV* 31 



24+0.6 
O.oM 



-0 4 +0 - 4 
>3.0 

i 1+0.8 
-0.7 



<10" 2 

4.0±1.5 
<l(r 2 



1 +20 



NOTE. — Col. (2): Galactic column density, calculated adopting the neutral hydrogen data set compiled by Dickey & Lockman (1990), using the CIAO observing toolkit at 
http://asc.harvard.edu/toolkit/colden.jsp. Col. (3): ACIS exposure time. Col. (4): total 0.5-10 keV counts. Col. (5): soft 0.5-2 keV counts. Col. (6): hard 2-10 keV counts. Col. (7): 
Hardness ratio, HR= (H -S)/(H + S). Col. (8): photon index of a power-law model where n(E) oc £™ r , assuming Galactic column density. Col. (9): Intrinsic column density estimated 
assuming a power4aw model with V = 1.8, atypical value for unobscured AGNs (e.g., Green et al. 2009), absorbed by a gas column /V H . 



ing to the inferred background counts level using the tables in 
Gehrels (1986) appropriate for small numbers of events. 

3.2.3. Source Decomposition for Marginally Resolved Nuclei 

The F-band nuclei of SDSS Jl 108+0659 and SDSS 
Jl 131-0204 are both separated by 0."70 (Table 1), which is 
close to the limit of the resolving power of ACIS. Neither one 
of the nuclei in SDSS Jl 131-0204 was detected in our obser- 
vations, and we estimated a 3-cr upper limit using the tables 
in Gehrels (1986) at the F-band position of each nucleus. 

To decompose the X-ray emission from the double nuclei in 
SDSS Jl 108+0659, we first performed a 2D image fitting in 
Sherpa, using PSF models as convolution kernels. We were 
unable to unambiguously separate the double nuclei, whose 
counts were too few for a statistically significant 2D decom- 
position. To increase the S/N, we then performed a ID anal- 
ysis for the nuclear X-ray emission in SDSS Jl 108+0659. 
There were too few counts to do a similar ID analysis in 
SDSS Jl 131-0204. 

For the ID analysis in SDSS Jl 108+0659, we projected the 
X-ray (full band) source profile in a direction connecting the 
two F-band nuclei by extracting the counts in a grid of re- 
gions of 0."5 x 2."0, as shown in Figure 2. The grid was 
designed to minimize contamination from the extended emis- 
sion to the northeast of the nuclear region. Figure 2 shows 
the projected ID source profile. We compared the ID pro- 
file with the 2D PSF model projected and convolved with the 
same grid. We performed fitting in Sherpa using x 2 statistics 
with the Gehrels variance function (Gehrels 1986); we tested 
two scenarios using one-component and two-component (i.e., 
two AGNs) PSF models, respectively. Centroids and ampli- 
tudes of the PSF models have been left free to vary. A likeli- 
hood ratio test suggests that the data favor a two-component 
marginally over a one-component model at the ~ 1 .5-a sig- 
nificance level. 

We show in Figure 2 our best-fit model for the ID nuclear 
source profile of SDSS Jl 108+0659. The best-fit central po- 
sitions of the two sources are consistent (within l-a uncer- 
tainties, 0."2) with the F-band positions projected on the axis 
connecting the two nuclei. We report the best-fit X-ray nu- 
clear positions of SDSS Jl 108+0659 in Table 2. We then re- 
peated the decomposition analysis for the soft and hard bands, 
respectively. With fewer counts in the individual band fittings, 
we fixed the centroids of the two models at the best fit posi- 
tions from the full band analysis, allowing only their ampli- 



tudes to vary. Figure 3 displays the best-fit models for the soft 
and hard band, respectively. The NW component is stronger 
(weaker) than the SE component in the soft (hard) band. The 
spectral properties of the two sources seem to be significantly 
different from each other, lending further support to the two- 
component scenario. In Table 4 we report the number of back- 
ground subtracted counts separately in the soft and hard bands 
as well as in the full band. 

3.2.4. X-ray Hardness Ratio and Spectral Analysis 

There are too few X-ray counts of our targets to perform 
reliable spectral fitting. Instead, we estimate parameters of X- 
ray spectral models using hardness ratios (HR) as a proxy for 
detailed spectral fitting. The hardness ratio is defined as 

where H and S are the number of counts in the hard and soft 
bands, respectively. We adopted the Bayesian Estimation of 
Hardness Ratios (BEHR; Park et al. 2006) to measure the HRs 
and their uncertainties, appropriate for the low count regime. 
To estimate the photon index T, we assume a single power- 
law model n(E) oc E~ r , absorbed by the Galactic column den- 
sity, calculated using the CIAO observing toolkit 14 based on 
the neutral hydrogen data set compiled by Dickey & Lock- 
man (1990). We also estimate the intrinsic galactic column 
density by fixing T = 1.8, typical for low-redshift unobscured 
Seyferts (Green et al. 2009). We list in Table 4 our X-ray 
count measurements and estimates for model parameters for 
each detected nucleus. 

We caution that the adopted single absorbed power-law 
model is most likely too simple for the X-ray spectra of ob- 
scured AGNs, in which thermal emission from starburst com- 
ponents and scattered nuclear emission are often present (e.g., 
Turner et al. 1997a,b). However, the low counts of our detec- 
tions do not allow us to test more realistic models. In addi- 
tion, our estimates of the intrinsic absorbing column may not 
necessarily reflect the true values in cases of patchy obscu- 
ration and/or significant scattering off an ionized medium in 
Compton-thick (i.e., /V H ~ 10 24 cm" 1 or larger) AGNs [which 
represent about half of the local Type 2 Seyfert population 
(Risaliti et al. 1999)], as observed in NGC 6240 (e.g., Vignati 
et al. 1999; Ptak et al. 2003) and in NGC 1068 (e.g., Matt et al. 

14 http://asc.harvard.edu/toolkit/colden.jsp. 
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FIG. 2. — Left panel: The grid region for constructing the ID spatial profile for SDSS Jl 108+0659. Background is the Chandra ACIS X-ray full-band image 
(unsmoothed with 0."25 pixel -1 binning). The blue circles indicate the Y-band positions and uncertainties of the double stellar nuclei from our HST WFC3 
imaging. Right panel: Projected ID profile of the nuclear X-ray emission in the full band. The grid positions correspond to those as shown by the grid region in 
the left panel and are separated by 0."5 (major tickmarks). Data are shown in diamonds with error bars, whereas our best-fit models are displayed in solid (for 
the total) and dotted/dashed (for each individual component) curves, respectively. The vertical lines indicate the X-ray centers of the two components. The blue 
circles on top indicate the projected 7-band positions of the double stellar nuclei. 




Same as in the right panel of Figure 2, but for the soft (left panel) and hard (right panel) band, respectively. 



FIG. 3.— 

1997; Guainazzi et al. 1999), although, again, the quality of 
our data do not allow us to robustly test these possibilities. 

4. RESULTS 

We examine the nature of the ionizing sources in the four 
optically selected kpc-scale binary AGNs. We first discuss 
the intrinsic X-ray luminosity and spectral properties (Section 
4.1). We then estimate the contribution from star formation 
related processes in the nuclear region to the observed X-ray 
flux (Section 4.2). Finally, we address whether these new X- 
ray and HST observations support the binary AGN scenario 
for each target (Section 4.3), estimate the X-ray-to- [Om] lu- 
minosity ratio of optical kpc binary AGNs, and compare with 
single AGNs (Section 4.4). 



4.1. X-ray Luminosity and Spectral Properties 

X-ray emission provides the most direct evidence for nu- 
clear activity. In particular, the 2-10 keV hard X-ray band is 
transparent to column densities of Nh < 10 24 cm -2 . To infer 
X-ray luminosity, we assume a simple absorbed power-law 
model, with a fixed T = 1.8 and the estimated intrinsic host 
galaxy column density as given in Table 4. In Table 5 we list 
the observed X-ray flux and luminosity as well as the esti- 
mated unabsorbed/intrinsic X-ray luminosity of each nucleus 
in the total, soft, and hard bands, respectively. 

The nuclei of our targets are optically classified as Type 2 
AGNs (Liu et al. 2010a), whose observed [Om] luminosities 
(Table 2) suggest moderate AGN luminosities. The four hard 
X-ray detected nuclei have estimated unabsorbed 2-10 keV 
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TABLE 5 
X-ray Fluxes and Luminosities 





A ,U . J — 1 UKC V 


Fv 1) ^ 9VpV 
A ,U. J — ZRcV 


Fv 9 i ni-pV 

A ,Z 1 UKC V 


Lv n ^ intpV 

,U.-J — 1 UK.C V 

dO 42 eres s~'l 


^\ ,U.J — ZK.C V 


Ly 9 intpV 

- B - J A ,Z — 1 UK.C V 


Ohiert Name 


("|Cy~ 14 pfpc c — 1 prn - ^ 


/in— 14 pi-nc c — 1 prn"2i 


no -14 eres s _1 cm" 2 ) 


HO 42 eres s" 1 ") 


HO 42 eres <T : 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SDSS J1108+0659NW 


1.05 ±0.15 


0.41±0.08 


0.46±0.23 


0.91±0.17 


0.34±0.06 


0.40±0.2() 


SDSS Jl 108+0659SE (observed) 


1.20±0.35 


0.054±0.030 


1.11±0.36 


1.08±0.31 


0.049±0.032 


1.00±0.29 


SDSS Jl 108+0659SE (unabsorbed) 


2.10±0.48 


0.59±0.38 


1.31 ±0.43 


1.89±0.53 


0.53±0.35 


1.18±0.42 


SDSS J1131-0204W 


<0.23 






<0.1 






SDSS J1131-0204E 


<0.23 






<0.1 






SDSS J1146+5110SW (observed) 


2.49±0.78 


0.08±0.07 


2.29±0.81 


1.05±0.33 


0.030±0.025 


0.96±0.34 


SDSS Jl 146+51 10SW (unabsorbed) 


5.2±1.6 


2.18±0.77 


3.2±1.1 


2.20±0.69 


0.92±0.28 


1.30±0.46 


SDSSJ1146+5110NE 


0.23±0.17 


0.23±0.17 


<0.9 


0.10±0.07 


0.10±0.07 


<0.30 


SDSS J1332+0606SW 


<0.35 






<0.37 






SDSS J1332+0606NE (observed) 


0.39±0.14 


0.062±0.050 


0.33±0.28 


0.43±0.16 


0.069±0.057 


0.37±0.30 


SDSS J1332+0606NE (unabsorbed) 


0.56±0.21 


0.21±0.16 


0.35±0.30 


0.65±0.24 


0.23±0.19 


0.39±0.31 



NOTE. — Cols. (2)-(4): total, soft, and hard X-ray flux and 1 a error. Col. (5)-(7): total, soft, and hard X-ray luminosity and 1 a error. For obscured sources, both observed (obs) 
and unabsorbed (unabs) estimates are listed. The unabsorbed estimates were calculated assuming a power-law model with V = 1.8 and A^h as derived from the HR measurements, as 
listed in Table 4. All errors quoted are statistical uncertainties only. 



luminosities ranging from 3.9 ±3.1 x 10 41 ergs s _1 to 1.3 ± 
0.5 x 10 42 ergs s" 1 , and unabsorbed 0.5-10 keV luminosities 
ranging from 6.5 ±2.4 x 10 41 ergs s" 1 to 2.2 ±0.7 x 10 42 ergs 
s -1 . The estimated upper limits for the four hard X-ray unde- 
tected nuclei range from ~ 1 .0 to 3.7 x 10 41 ergs s -1 in 0.5-10 
keV. These luminosity estimates are similar to the few previ- 
ously known X-ray confirmed kpc binary AGNs [NGC 6240 
(Komossa et al. 2003), 3C 75 (Hudson et al. 2006), Mrk 463 
(Bianchi et al. 2008), Mrk 266 (Brassington et al. 2007), and 
Mrk 739 (Koss et al. 2011); see also Ballo et al. (2004) for 
a candidate Arp 299]. We devote the rest of the section to 
determining the nature of the X-ray sources. 

First, we consider the possibility of ultra-luminous X-ray 
sources (ULXs, with typical X-ray luminosities of order 10 39 - 
10 40 ergs s" 1 , i.e., beyond high mass X-ray binaries but much 
less than typical AGNs; Long & van Speybroeck 1983; Fab- 
biano 2006), which are off-nuclear point-like X-ray sources 
commonly observed in local major mergers of disk galaxies 
such as in Arp 244 (i.e., the Antennae) and Arp 270 (e.g., 
Brassington et al. 2007). Some ULXs may be accreting in- 
termediate mass BHs (M BH ~ 10 2 -10 4 M Q ; e.g., Miller et al. 
2003; Miller & Colbert 2004). Since the nuclei in our tar- 
gets have X-ray luminosities much higher than typical ULXs, 
they are most likely of a different origin. Further evidence 
against the ULX scenario includes: (1) the BHs are expected 
to be supermassive given their host bulge properties (Liu et al. 
2010a); (2) the observed [Om] luminosities are significantly 
higher than those for ULXs (Abolmasov et al. 2007); and (3) 
the X-ray point sources in our targets are nuclear given astro- 
metric uncertainties, although ULXs could also live close to 
nuclear regions. 

Second, the estimated intrinsic hard X-ray luminosities of 
our targets are close to or below ~ 10 42 ergs s" 1 - the charac- 
teristic upper limit for the most luminous star-forming galax- 
ies (e.g., Zezas et al. 2001). So it is quite possible that much 
or all of the luminosity is due to star formation. While X-ray 
spectral shape offers another diagnostic to discriminate be- 
tween AGN and starburst scenarios, the uncertainties of our 
spectral estimates are too large to draw firm conclusions for 
the majority of the nuclei. Therefore we must factor in some 
independent star formation rate estimates to critically test the 
AGN scenario for each nucleus. 



4.2. Contribution from Nuclear Star Formation 

Even without the presence of an AGN, star formation- 
related processes [e.g., accretion onto a neutron star or a 
black hole in X-ray binaries, thermal bremsstrahlung from a 
starburst-driven wind (e.g., Strickland & Stevens 2000)] may 
produce strong soft and hard X-ray emission. Intense nuclear 
star formation often accompanies and sometimes outshines 
AGN in X-rays in gas-rich mergers, making kpc-scale binary 
AGNs challenging to pin down. To break the degeneracy, we 
first examine the contribution from nuclear star formation, us- 
ing independent constraints from HST U -band imaging. 

We estimate the expected X-ray emission due to star for- 
mation within the same apertures used to perform our X-ray 
extraction, to evaluate if an additional excitation source, i.e., 
an AGN component, is needed. To derive X-ray luminosites 
from SFRs, we adopt the empirical calibration of Ranalli et al. 
(2003, see also Grimm et al. 2003) based on 23 nearby star- 
forming galaxies, which is given by 

Af 5 -2kev = 4.5 x 10 39 -^ergs s" 1 , (3) 
Mq yr 

il-iokev = 5-0 x 10 39 -^ rr ergs s" 1 , (4) 
m & yr 

with an rms scatter of 0.27 dex and 0.29 dex, respectively. In 

Table 3 we list the derived ^o I 5_2keV ^2-iokeV estimates for 
each nucleus. The predicted X-ray contribution from star for- 
mation is an order of magnitude or more below the observed 
X-ray luminosity. At this point, the case that the X-rays are 
coming from AGN seems unambiguous, before considering 
the uncertainties. 

4.2.1. Uncertainties 

In this section, we discuss systematics and uncertainties of 
our estimates of the expected X-ray luminosity due to star 
formation-related processes. First, a Salpeter initial mass 
function (IMF; Salpeter 1955) was assumed with mass limits 
of 0. 1 Mq and 1 00 M in the adopted calibration of SFR from 
f/-band luminosity. The adopted SFR-X-Ray-luminosity re- 
lation of Ranalli et al. (2003) was calibrated under the same 
assumptions about the IMF and mass range (essentially all 
from Kennicutt 1998), so that the inferred SF-related X- 
ray luminosity is not sensitive to the IMF uncertainty, given 
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FIG. 4. — X-ray luminosities versus the expected contribution from star formation-related processes. We show the soft-band results in the left panel and the 
hard-band results in the right panel, respectively. Open symbols represent observed quantities, whereas filled symbols are corrected for dust extinction and 
gas absorption. We use dotted lines to connect values of the same nucleus before and after extinction/absorption correction. For comparison we also show 
measurements of NGC 6240 (total emission from both nuclei) from the literature (Ptak et al. 2003). In both panels, we show the equality relation with a solid 
line. Error bar in the lower right of the panel indicates typical uncertainty (~ 0.4 dex) in the X-ray luminosity from star formation-related processes, which was 
estimated by convolving that propagated from SFR estimates with the rms scatter of the SFR-Lx calibrations. 

in the soft and hard bands, respectively. Armed with both 
the observed X-ray properties and independent constraints on 
the expected X-ray contribution due to star formation from 
[/-band imaging, we now discuss whether the observations 
support the binary-AGN scenario for each of our targets. 

4.3.1. SDSS J 1108+0659 

Both nuclei were optically classified as Type 2 Seyferts, 
suggesting the presence of at least one AGN component. Both 
nuclei were detected in both soft and hard X-ray bands. The 
[/-band image reveals intense star formation activity in the 
NW nuclear region. While the X-ray HR and T estimates 
(based on the simple absorbed power-law model) suggest no 
nuclear obscuration, the strong starburst component may be 
accompanied by significant dust and gas, suggesting the pres- 
ence of a substantial absorbing column. The adopted sin- 
gle absorbed power-law model is most likely too simple, but 
there are too few X-ray counts to test more realistic multi- 
component spectral models. The spectral properties of the 
SE source suggest moderate nuclear obscuration, with an es- 



that there is no large systematic IMF variations among star- 
forming galaxies (Scalo 1986; Kroupa 2001; Chabrier 2003). 

The uncertainty of the estimated X-ray luminosity due to 
star formation is likely to be dominated by the poorly con- 
strained [/-band extinction correction. First, we estimated 
Fho/Fhp using emission line measurements for each nucleus 
from our ground-based optical longslit spectroscopy. Due to 
its lower angular resolution and projection effects, the aper- 
ture of the emission line measurement does not exactly match 
with the [/-band measurement. To estimate the reddening 
uncertainty due to this coverage mismatch, we have com- 
pared the emission-line ratio measurement from our slit spec- 
troscopy against that from the SDSS fiber spectra. For the 
SDSS fiber-integrated measurements, we assumed that the 
two velocity components correspond to the two nucleus com- 
ponents, although the actual association is likely to be more 
complicated. Nevertheless, we found that the uncertainty in 
[/-band extinction estimates due to aperture mismatch is in 
general < 1 mag (Table 3). Multi-color imaging and/or IFU 
spectroscopy with the same angular resolution as the [/-band 
imaging would help further constrain the uncertainty. 

More importantly, the color excesses we derived using the 
Balmer decrement method do not necessarily represent the 
true dust content of our targets. For example, most of the 
dust could be concentrated on scales smaller than where the 
Balmer lines are emitted. High extinction with little redden- 
ing could arise in objects with very patchy and optically thick 
dust clouds. These caveats associated with the uncertainties 
in the extinction geometry and the reddening assumptions are 
typical, but must be kept in mind when extinction corrections 
are applied to compare with other studies in the literature. 

4.3. Nature of the Ionizing Sources 

In Figure 4 we compare the expected X-ray luminosities 
due to star formation against the observed X-ray luminosities 



timated column density Nh ~ x 10 cm . The SE nucleus 
was undetected in the U -band, indicating very low level of star 
formation around the SE nucleus (or that the star formation 
was highly obscured). The spatial profiles of the two nuclear 
X-ray sources support the AGN scenarios for both, although 
the starburst components in the NW nucleus are too compact 
to be resolved in the X-rays. For both nuclei, the expected SF- 
related X-ray luminosities are too low (5.7 ± 1 .0 times fainter 
in the soft and 6.7 ± 3.3 times fainter in hard X-ray band for 
the NW nucleus, and at least three orders of magnitude fainter 
in both bands for the SE nucleus) to explain the observed val- 
ues, lending further support to the binary-AGN scenario. 

4.3.2. SDSS Jl 131 -0204 
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Both nuclei were optically classified as Type 2 Seyferts, in- 
dicating that at least one AGN component must be present. 
Both nuclei were undetected in the two X-ray bands. The 
merger system seems to be embedded in a massive disk com- 
ponent, and we may be viewing the disk at an angle close to 
edge-on through a large amount of absorbing column. Our U - 
band image reveals circumnuclear star formation around the 
eastern nucleus, whereas the western nucleus was undetected 
in the J/-band. The U -band constraints on the expected SF- 
related X-ray luminosities are still roughly an order of mag- 
nitude smaller than the observed X-ray luminosity upper lim- 
its for both nuclei. Therefore, the upper limits of the X-ray 
luminosities are still consistent with the presence of double 
AGNs, although the possibility that one AGN ionizes gas in 
both galaxies cannot be ruled out. 

4.3.3. SDSSJ1146+5110 

The SW nucleus was optically classified as a Type 2 
Seyfert, whereas the NE nucleus is either a Type 2 Seyfert, 
a LINER, or a LINER-H h composite. The X-ray spatial pro- 
files of the two nuclei support the AGN scenario for both, 
although both nuclear starburst components are too compact 
to be resolved in the X-rays. The SW nucleus was detected 
in both soft and hard X-rays. Its HR and T measurements 
may suggest a mild obscuration (Nh ~ 4.0 ± 1.5 x 10 22 cm" 2 
estimated assuming T=1.8); its intrinsic hard X-ray luminos- 
ity (Lx ; 2-iokev ~ 1.3 ±0.5 x 10 42 erg s" 1 ) is over two or- 
ders of magnitudes larger than that expected from star for- 
mation (Figure 4), strongly suggesting an AGN component. 
The NE nucleus was detected only in the soft X-rays. Taken 
at face value, its HR and T estimates suggest a steep spec- 
trum; this, together with its moderate hard X-ray luminos- 
ity (iz,2-iokeV < 3.0 x 10 41 ergs s" 1 ) may indicate a source 
dominated by star formation-related processes. However, the 
expected star formation-related X-ray luminosities are lower 
than the observed value (14 ± 10 times fainter in the soft 
band) or the upper limit (up to ~ 40 times fainter in the hard 
band), lending support for an AGN component also in the 
NE nucleus. Similar to the case of the NW nucleus of SDSS 
J 1 108+0659, the U -band image reveals intense star formation 
activity in the NE nuclear region of SDSS Jl 146+5 1 10. While 
the apparent X-ray HR and T estimates suggest no significant 
absorbing column, the nuclear starburst component may indi- 
cate otherwise. 

4.3.4. SDSS J 1332+0606 

Both nuclei were optically classified as Type 2 Seyferts. 
The NE nucleus was detected in both soft and hard X-rays, 
although the counts were very low. The apparent HR and V 
measurements may suggest a reflection dominated spectrum, 
indicating the presence of a Compton-thick source. In the 
Compton-thick scenario, the true intrinsic X-ray luminosities 
would be much higher than our fiducial estimates, where only 
modest absorption is assumed (A^h ~ 1 -O^' 99 x 10 22 cm" 2 esti- 
mated assuming T=1.8, i.e., consistent with zero absorption). 
The expected SF-related X-ray luminosities in both bands are 
lower than the observed values within the uncertainties, sug- 
gesting an AGN component. The SW nucleus was undetected 
in the X-rays. While the upper limits of the X-ray luminosi- 
ties are consistent with an AGN component in the S W nucleus 
as well, the possibility that only one AGN in the NE nucleus 
ionizes gas in both galaxies cannot be ruled out. 

In summary, our new Chandra and HST observations sup- 



port the binary-AGN scenario for two of our four targets 
(SDSS Jl 108+0659 and SDSS Jl 146+51 10). For the other 
two targets (SDSS Jl 131-0204 and SDSS J 1332+0606), the 
existing data are still consistent with the binary-AGN sce- 
nario, although the possibility of only one AGN ionizing both 
components in the mergers cannot be ruled out. 

4.4. X-ray to [O in] Luminosity Ratio 

For optically selected Type 2 AGNs, the [0 111] emission- 
line luminosity is usually taken as a surrogate to estimate the 
intrinsic hard X-ray luminosity, because the intrinsic AGN 
continuum luminosity is obscured in the optical. Measure- 
ments of the ratio Lx,2-iokev /£[Oiii] f° r optically selected Type 
2 AGNs span a wide range, with values ranging from a few to 
a few hundred (Mulchaey et al. 1994; Heckman et al. 2005; 
Panessa et al. 2006). 

4.4. 1 . Systematically Smaller X-ray-to-[ O III] -luminosity Ratio in 
Optically Selected kpc-Scale Binaries than in Single AGNs 

In Figure 5 we plot the hard X-ray luminosity against the 
[Oiii] luminosity for each nucleus in our targets. For con- 
text, we compare our targets to observations of single AGNs 
and the few previously known kpc-scale binary AGNs. We 
study both the relation between the observed hard X-ray lu- 
minosity Lx.2-iokeV,observed and the observed [Oiii] luminos- 
ity L[o in], observed, and that between the unabsorbed hard X-ray 
luminosity Lj^-iokev.unabsorbed and the extinction-corrected 
[Oiii] luminosity L[oni] extinction-corrected- We used different 
comparison samples for the two as usually only the observed 
or the corrected luminosity was available in any given sample. 

For the LA-^-iokeV^bserved-^iomi.observed relation (left panel of 
Figure 5), the comparison samples include the 47 hard X-ray 
(3-20 keV, in this particular case) selected AGNs (the z < 0.2 
subset of the Sazonov & Revnivtsev 2004 sample from the 
RXTE all sky survey, Revnivtsev et al. 2004) and 55 optically 
selected local [Om]-bright AGNs (Xu et al. 1999; Whittle 
1992) studied by Heckman et al. (2005), and 8 optically se- 
lected Type 2 quasars from Ptak et al. (2006). The Heckman 
et al. (2005) local AGN samples have similar redshifts and 
£[Oiii],observed to our targets, whereas the Ptak et al. (2006) ob- 
jects are at higher redshifts (z ~0.3-0.8 compared to our tar- 
gets at z ~0.1-0.2) and have higher L[om], observed (by ~ 1.5 
dex). Heckman et al. (2005) showed that optically selected 
Type 2 AGNs have systematically lower L;f,2-iokev,observed (by 
an average of 1.0 dex) at a given L[onrj, observed than hard X- 
ray selected AGNs (both Type 1 and Type 2) and optically 
selected Type 1 AGNs. Our optically selected kpc-scale bi- 
nary AGNs seem to have systematically smaller Lx/^[0iii] 
(observed) values than even optically selected single Type 2 
AGNs, although the sample size is still too small to draw a 
firm conclusion. The observed hard X-ray luminosities of our 
targets on average are <~ 0.8 ± 0.2 dex smaller at least 15 than 
those of optically selected single-nucleus Type 2 AGNs by 
Heckman et al. (2005) at a fixed Z40111], observed (i-e-, average 
log(Lx/£[Oiii]) = -0.2 with an rms scatter of 0.40 dex, com- 
pared to \og(L x / L[om]) = 0.6 with an rms scatter of 1.1 dex 
for optically selected single Type 2s). 

For the Lx,2-10keV,unabsorbed-£[Om]extincuon-corrected relation 

(right panel of Figure 5), the comparison samples include 47 
Palomar Seyfert galaxies (optically selected Type 1 and Type 
2 Seyferts drawn from the Palomar survey of nearby galaxies 

15 The upper limits were included in the fit as detections. 
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FIG. 5. — Hard X-ray versus [OIII] A5007 emission line luminosity. Our targets are indicated with large colored open symbols, as denoted in the plot. Left 
panel: observed 2-10 keV luminosity versus observed [OIII] luminosity. The comparison samples include: hard X-ray selected AGNs (Type Is as filled and 
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green dash-dot-dot line is for the CSC-SDSS sample, whereas the cyan dotted line is that from Panessa et al. (2006) for mixed Seyferts in nearby galaxies. Shown 
as filled bowties are two NLR-kinematics-candidate double-peaked [OIII] AGNs [left panel: SDSS J1715+6008 from Comerford et al. (2011), and right panel: 
CXO J1426+35 (total emission) from Barrows et al. (2012)]. 



by Ho et al. 1995) from Panessa et al. (2006), and three pre- 
viously known kpc-scale binary AGNs [NGC 6240 from Ko- 
mossa et al. (2003); Ptak et al. (2003), Mrk 463 from Bianchi 
et al. (2008), and Mrk 266 from Brassington et al. (2007)]. 
We also include a new sample of 55 single-nucleus Type 2 
Seyferts from Trichas et al. (2012). This comparison sample 
was both optically and X-ray selected, and was constructed 
by cross-matching the Chandra Source Catalog (CSC; Evans 
et al. 2010) Release 1.1 with the SDSS DR7 spectroscopic 
galaxy catalog at z < 0.3, and selecting objects whose opti- 
cal emission line ratios [Om] A5007/H/? and [Nil] A6584/Ha 
are characteristic of Type 2 Seyferts according to the Kewley 
et al. (2001) criteria. 

Panessa et al. (2006) suggested that after properly account- 
ing for absorption correction (including for Compton thick 
sources), optically selected Type 1 and Type 2 Seyferts all 

follow the Same Lx,2-10keV,unabsorbed-i[Oni].extinction-coiTected re- 
lation. In particular, optically selected Type 2 Seyferts, which 
were significantly X-ray weaker than Type 1 Seyferts and 
quasars, also obey the same relation, after the "Compton 
thick" luminosity correction. The average absorption correc- 
tion for Compton-thick sources (30% of the sample) is ^2-3 
dex, or ^0.5-1 dex for the other Type 2 Seyferts. Their ab- 
sorption correction, however, is significantly larger than that 
for our targets, although our correction may have been under- 
estimated given the poor X-ray constraints. 

After correction for gas absorption and dust extinction, the 
unabsorbed hard X-ray luminosities of our targets appear to 
be ~ 2.4 ±0.3 dex smaller (at logLjomj of 42.0) than those 



expected from the Panessa et al. (2006) relation, logLx = 
1 .221ogL[oni]-7.34, although the absorption correction of our 
targets may have been significantly underestimated (Section 
3.2.4). Our targets have average unabsorbed log(Lx/i[om]) = 
-0.5 (with an rms scatter of 0.69 dex), which is ~ 1.9 ±0.3 
dex lower than that of the CSC-SDSS cross-match com- 
parison sample (average log(Lx/i[Oiii]) =1-4 with an rms 
scatter of 0.77 dex). Like our targets, the three previously 
known kpc-scale binary AGNs (all of which are ULIRGs and 
were discovered in the X-rays) also show smaller unabsorbed 
log(Lx/L[om]) values than the Panessa et al. (2006) relation. 

4.4.2. Interpretation: Higher Nuclear Gas Concentration and/or 
Viewing Angle Effect 

We now discuss possible causes of the apparent hard X-ray 
weak tendency in optically selected kpc-scale binary AGNs 
compared to optically selected single Type 2 AGNs. Sim- 
ulations suggest that galaxy mergers may funnel significant 
amount of gas towards galaxy centers, triggering both global 
and nuclear starburst activity as well as AGN (e.g., Hernquist 
1989; Hopkins et al. 2008). Using a sample of 1286 AGN 
pairs spectroscopically selected from the SDSS with projected 
separations of a few kpc to a few tens of kpc (Liu et al. 201 1), 
Liu et al. (2012b) have shown that the fraction of both single 
and double AGNs stays constant with decreasing projected 
separation at > 20 kpc scales. This suggests that at wide sep- 
arations, the majority of the observed AGNs are most likely 
due to stochastic accretion not associated with tidal interac- 
tions (e.g., Ciotti & Ostriker 2007; Ciotti et al. 2010). How- 
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ever, the fraction of both single and binary AGNs increases 
rapidly on scales below 20 kpc, indicating tidally enhanced 
AGN in close galaxy pairs (see also Ellison et al. 2011). In 
addition, the fraction of binary AGNs in close galaxy pairs is 
significantly higher than what would be expected from ran- 
dom pairing of single AGNs. This, together with the correla- 
tion observed in AGN luminosity between the merging com- 
ponents in a pair (Liu et al. 2012b), provides direct evidence 
that the double AGN observed in close galaxy pairs may in- 
deed be connected to merger activity. It is likely that the gas 
in kpc-scale binary AGNs is being tidally funneled all the way 
to the nuclear regions of both components, and either fuels the 
binary AGN directly, or fuels a nuclear starburst which in turn 
feeds the BHs. 

We therefore suggest that a higher gas column in the nuclear 
regions of kpc-scale binary AGNs, which is likely induced by 
the merger events, may be at the root of the observed X-ray 
weak tendency. Unlike X-ray emission which comes from the 
accretion disk/corona of the accreting SMBHs, [O in] emis- 
sion in the NLRs comes from much larger scales (Antonucci 
1993; Urry & Padovani 1995), and should therefore be much 
less subject to nuclear obscuration. In this scenario, the X- 
ray absorption column was significantly underestimated for 
our targets, which is plausible given the large systematic un- 
certainty of our X-ray measurement. While the three previ- 
ously known kpc-scale binary AGNs (which are all ULIRGs 
and whose binary- AGN nature was all discovered in the X- 
rays) also show smaller X-ray-to- [Om] luminosity ratios than 
single Seyferts, the X-ray weak tendency seems to be less 
pronounced than for our optically selected kpc-scale binary 
AGNs. This may be due either to a selection effect, or a host- 
galaxy effect (i.e., the nuclear gas distribution in ULIRGs - 
disk-dominated mergers - is less concentrated than in opti- 
cally selected binary AGNs - whose hosts have prominent 
stellar bulges - due to lower bulge-to-disk ratios of the host 
galaxies), or a combination of both. 

Alternatively, the X-ray weak tendency observed in our tar- 
gets may be caused by a viewing angle effect related to its 
double-peak selection. The requirement of detecting well- 
separated [Om] double peaks in velocity is likely to pick 
out edge-on systems, and therefore would be biased towards 
a higher absorbing gas column. For example, SDSS J 1 131- 
0204, the only one of our targets which was undetected in 
the X-rays, appears to show two stellar nuclei imbedded in a 
large edge-on disk (Liu et al. 2010a; Shen et al. 2011). To 
check whether the double-peak selection causes the X-rays to 
be weak, we also examine two double-peaked [O in] AGNs in 
the literature with available [O in] and hard X-ray luminosity 
measurements [SDSS J1715+6008 from Liu et al. (2010b), 
Smith et al. (2010), and Comerford et al. (2011), and CXO 
J1426+35 from Barrows et al. (2012)], whose emission line 
profiles are more likely caused by complex NLR kinemat- 
ics, i.e., not caused by kpc-scale binary AGNs. As shown in 
Figure 5, in both cases, the X-ray-to- [Om] ratios are smaller 
than the Heckman et al. (2005) and Panessa et al. (2006) re- 
lations for the observed and corrected luminosities, respec- 
tively. While this comparison is based on only two objects, it 
may suggest that the hard X-ray weak tendency is related to 
a double-peak selection, regardless of whether there are two 
active black holes in the system. However, the viewing angle 
effect is unlikely to be the sole cause of the X-ray weak ten- 
dency, given that the three previously known kpc-scale binary 
AGNs, which were not selected from double-peaked narrow 
line AGNs, also show smaller X-ray-to- [Om] luminosity ra- 



tios than single Seyferts. 

In summary, we conclude that the observed X-ray weak 
tendency in our optically selected kpc-scale binary AGNs is 
likely caused by a combination of: (a) a higher nuclear gas 
column, which may be induced by merger events, and (b) a 
viewing angle bias related to its double-peak narrow emission 
line selection. Of course, another possibility is that we were 
simply unlucky with these four systems, and they are not typ- 
ical of kpc-scale binary AGNs. 

5. DISCUSSION 

We discuss the implications of our results for the gen- 
eral double-peak narrow line selection approach of identify- 
ing kpc-scale binary AGNs (Section 5.1), for the attributes 
and limitations of optical identification compared to X-ray 
searches (Section 5.2), and for the observed frequency of kpc- 
scale binary AGNs (5.3). 

5.1. On the Double-Peak Approach for kpc-Scale 
Binary -AGN Identification 

5.1.1. X-ray Confirmation of Optical Candidates: Success and 
Ambiguity 

Until recently, searches for kpc-scale binary AGNs have 
been more or less serendipitous. Previous work (Liu et al. 
2010a; McGurk et al. 2011; Shen et al. 2011; Fu et al. 
2012) has demonstrated that selecting spectroscopic candi- 
dates based on double-peaked narrow emission lines and 
follow-up with high spatial resolution imaging (Fu et al. 
2011a; Rosario et al. 2011) and spatially resolved spec- 
troscopy (Comerford et al. 2012) is a promising technique to 
identify genuine kpc-scale binary AGNs. With the new Chan- 
dra and HST observations presented here, we have critically 
examined the nature of the ionizing sources of the optically 
identified kpc-scale binary AGNs. By combining SFR con- 
straints from high resolution U -band imaging with spatially 
resolved X-ray spectroscopic imaging, we put much stronger 
and more direct constraints on the ionization mechanism of 
the close double nuclei than using X-ray observations or op- 
tical diagnostics alone. Our results confirm the kpc-scale bi- 
nary AGN nature for two of the four optical candidates; the 
data are still consistent with the binary scenario for the other 
two, but we cannot rule out the possibility of one AGN ioniz- 
ing gas in both merging components. While the result lends 
some further support to the overall approach of systematically 
identifying kpc-scale binary AGNs based on the double-peak 
selection, it also suggests that X-ray confirmation of optical 
binary candidates can be challenging and ambiguous, at least 
for heavily absorbed sources. 

5.1.2. Importance of Identifying Double Nuclei Both in Gas and in 

Stars 

Using follow up observations, Shen et al. (2011, see also 
Fu et al. 2012) suggested that the majority of optically se- 
lected AGNs with double-peaked [Om] AA4959,5007 emis- 
sion lines are not due to kpc-scale binary AGNs, but are 
caused by complex NLR gas kinematics around single AGNs, 
such as outflows (e.g., Rosario et al. 2010) and/or rotating 
disks (e.g., Smith et al. 2012; Blecha et al. 2012b). Shen et al. 
(2011) have emphasized the importance of combining high 
resolution NIR imaging (to resolve the double stellar bulges 
expected to be associated with the BHs, and to detect tidal 
features which are commonly observed in major mergers with 
kpc-scale nuclear separations) with spatially resolved optical 
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spectroscopy (to locate the ionizing sources and to register 
with the BHs/bulges) to identify strong kpc binary AGN can- 
didates (see also McGurk et al. 2011). Detecting two stellar 
components in a galaxy merger alone (e.g., Fu et al. 2011a; 
Rosario et al. 201 1) is not a sufficient condition for identifying 
binary AGNs, as the double-peaked profile could be caused 
by NLR gas kinematics around a single AGN in one of the 
merger components (Shen et al. 2011). Neither is detecting 
spatial offsets between the emission line components alone 
(e.g., Gerke et al. 2007; Comerford et al. 2009a; Xu & Ko- 
mossa 2009; Barrows et al. 2012) a sufficient condition, even 
though most double-peaked AGNs show spatial offsets on at 
least kpc scales between the redshifted and blueshifted emis- 
sion line components (Shen et al. 201 1; Fu et al. 2012; Comer- 
ford et al. 2012), because such spatial offsets are expected and 
also commonly observed in single AGNs due to the spatial ex- 
tent of the NLR, as demonstrated by Fischer et al. (2011) in 
Mrk 78 (see also Whittle & Wilson 2004). 

To pin down the ionizing sources, candidates identified by 
spatially resolved optical spectroscopy can be confirmed us- 
ing imaging spectroscopy in the X-rays (as in this work) or 
in the radio (Fu et al. 2011b). However, identifying the dou- 
ble stellar bulges expected to be associated with the purported 
double BHs should be a prerequisite before carrying out X- 
ray nor radio observations. We emphasize the importance 
of detecting double nuclei in both gas and in stars in the 
identification of kpc-scale binary AGNs. X-ray and/or ra- 
dio observations would help clarify the ambiguities for gas 
(i.e., in [Oiii]), but cannot substitute deep NIR imaging for 
stars. Double X-ray and/or radio sources (as well as spatially 
resolved optical emission lines coincident with the two X- 
ray/radio sources) on a smooth stellar background (e.g., see 
the candidate reported by Comerford et al. 2011), i.e., with- 
out the expected double stellar bulges, is most likely evidence 
for complex NLR gas kinematics (e.g., a jet; Comerford et al. 
2011) rather than binary AGNs, as the case of Mrk 78. Ex- 
ceptions to this included binary AGNs with separations ^100 
pc, where the two stellar bulges may have merged, and minor 
mergers, which are difficult to resolve due to to large con- 
trast ratios (see the candidate example discovered by Fabbiano 
etal. 2011). 

5.1.3. Importance of High-Resolution Deep Observations 

Our new Chandra and HST observations here have unam- 
biguously confirmed the kpc-scale binary AGN nature for 
SDSS Jl 108+0659 and SDSS Jl 146+51 10. Fu et al. (2012) 
suggested that the case for SDSS Jl 108+0659 was ambigu- 
ous from the detection of extended emission line regions 
(EELRs), using Keck NIR imaging assisted with laser guide 
star AO and seeing-limited integral field spectroscopy (IFS) 
on the University of Hawaii 2.2-m telescope. However, the 
angular resolution of the IFS data (seeing ~ 0."8) may have 
been insufficient to fully resolve the close double nuclei in 
SDSS Jl 108+0659 (angular separation of 0."70 as measured 
by our HST imaging). More importantly, the detection of 
EELRs on larger scales should not be taken as evidence for or 
against the binary AGN scenario, because EELRs have been 
observed in both kpc-scale binary (such as in NGC 6240 and 
in Mrk 266) and single AGNs. 

Fu et al. (2012) also suggested that SDSS Jl 146+51 10 was 
a pair of EELRs powered by a single AGN, but their IFS data 
were not sensitive enough to detect the faint [O in] emission 
associated with the NE nucleus, which was detected by Liu 
et al. (2010a) using longslit spectroscopy with the APO 3.5-m 



telescope. 

5.2. Comparing Optical to X-ray Identification of kpc-Scale 
Binary AGNs 

We discuss our results in the context of comparing optical 
and X-ray identification of kpc-scale binary AGNs. A com- 
plete census of the AGN population needs multi-wavelength 
probes. The identification of AGNs based on any particular 
wavelength window is likely to be limited by biases and in- 
completeness. For example, Heckman et al. (2005) suggested 
that identifying AGN based on bright [Oiii] emission lines 
will uncover the majority of hard-X-ray selected AGNs 16 , 
whereas identifying AGN by hard X-rays will miss a signifi- 
cant population of AGNs selected based on the [O in] optical 
emission line. These hard-X-ray faint yet [O in] bright Type 
2 AGNs are generally interpreted as heavily obscured or even 
Compton-thick sources (e.g., Risaliti et al. 1999; Bassani et al. 
1999; Levenson et al. 2002). 

Hard X-ray searches for kpc-scale binary AGNs are still 
limited by low angular resolution and small number statistics. 
Low resolution studies [by Jimenez-Bailon et al. 2007 using 
the XMM-Newton (Jansen et al. 2001) and by Koss et al. 2012 
using the Swift Burst Alert Telescope (Barthelmy et al. 2005) 
All-Sky Hard X-ray Survey (Ajello et al. 2008; Tueller et al. 
2010)] are restricted to wide (tens-of-kpc scales) pairs in the 
local universe, whereas high resolution searches with Chan- 
dra are confined to small samples (Teng et al. 2005; Koss et al. 
2012; Teng et al. 2012). In addition, X-ray searches often 
have to apply some pre-selection based on optical AGN diag- 
nostics to boost the success rate (Koss et al. 2012; Teng et al. 
2012). We use the result of this paper and previously discov- 
ered X-ray identified cases to understand the relative merits 
of optical and X-ray identification of kpc-scale binary AGNs. 

As with heavily obscured single AGNs, our result on the 
hard-X-ray-to-[Oin] luminosity ratio suggests that identify- 
ing kpc binary AGNs by hard X-rays is likely to miss a pop- 
ulation of optically selected Type 2 binary AGNs. To address 
the fraction of such heavily X-ray-obscured binary AGNs 
among all binary AGNs would require knowledge of the space 
densities of Type- 1 -Type- 1 and Type- l-Type-2 binaries, both 
optically selected and hard X-ray selected, which is beyond 
the scope of this paper. But given our result that the hard- 
X-ray-to-[Om] A5007-luminosity ratio of optically selected 
Type 2 binary AGNs appears to be systematically smaller than 
single Type 2 AGNs, we speculate that the obscured fraction 
could even be higher in binary than in single AGNs. 

On the other hand, the parent sample of double-peaked 
narrow emission line AGNs (from which we identified the 
kpc-scale binary AGN candidates) was selected to be Type 
2 Seyferts (Liu et al. 2010b) according to the Kewley et al. 
(2001) criteria. This, by construction, will miss Type 1 binary 
AGNs 17 , mixed Type 1/Type 2, as well as Type 2 binaries in- 
volving LINERS and/or AGN-H h composites. In particular, 
the two nuclei of the prototypical kpc binary AGN NGC 6240 

16 Except for BL Lac objects (Wolfe 1978) which contain very weak opti- 
cal emission lines (Collinge et al. 2005). 

Our parent AGN sample was selected from the SDSS DR7 main galaxy 
sample (Strauss et al. 2002), supplemented with Type 2 quasars from Reyes 
et al. (2008) that were not included in the main galaxy sample. The main 
galaxy sample includes objects with redshifts z < 0.7 and spectral classifica- 
tion as galaxies by the specBS pipeline (Adelman-McCarthy et al. 2008), or 
quasars that were targeted as galaxies. While Hao et al. (2005) have shown 
that the narrow emission line ratios of Type 1 AGNs also follow the Kewley 
et al. (2001) criteria, our parent AGN sample will miss (luminous) Type 1 
AGNs which were not included in the DR7 main galaxies. 
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(Komossa et al. 2003) are both heavily obscured in the optical, 
with emission line ratios characteristic of LINERs (Lutz et al. 
1999); another example is Mrk 266, whose northern nucleus 
is optically classified as a composite yet does contain an X-ray 
AGN (Brassington et al. 2007; Mazzarella et al. 2012). Liu 
et al. (2011) have addressed the frequency of binary AGNs 
(including both Type 1 and Type 2 objects as well as LIN- 
ERs and composites) on a few kpc to tens-of-kpc scales using 
spectroscopic galaxy pairs from the SDSS DR7. Complemen- 
tary to the double-peak approach, the Liu et al. (201 1) sample 
includes both Type 1 and Type 2 AGNs but is more focused on 
tens-of-kpc scales and is incomplete for binaries closer than 
5 kpc in projection. While more work is clearly needed to 
address the frequency of kpc-scale Type 1 binaries and Type 
2 binaries containing LINERS/composites, the frequency of 
kpc-scale binary AGNs inferred from the Liu et al. (2011) 
sample is similar to that from the double-peak approach (Liu 
et al. 2010a; Shen et al. 2011), suggesting that the majority 
of kpc-scale binary AGNs are Type 2 objects, most of which 
are Seyferts at least in the moderate AGN-luminosity regime 
being considered. 

5.3. Frequency of kpc-Scale Binary AGNs 

Motivated by the serendipitous success in NGC 6240, Teng 
et al. (2005) carried out a systematic search for binary AGNs 
using Chandra imaging of a sample of eight (U)LIRGs which 
contain double stellar nuclei; they found no additional strong 
candidate. Similarly, Teng et al. (2012) carried out a Chan- 
dra surveys of a sample of 12 massive galaxy mergers, each 
of which contains one optical AGN, but found no convincing 
case of binary AGNs. The null results from such "blind" X- 
ray searches in small samples are not unexpected, considering 
the relatively low frequency of kpc binary AGNs as measured 
using optical identification based on much larger parent sam- 
ples (Liu et al. 2010a; Shen et al. 2011; Fu et al. 2012; Liu 
et al. 2011). 

Using follow-up NIR imaging and optical slit spectroscopy 
of 43 double-peaked [Om] Type 2 AGNs from the parent 
sample presented in Liu et al. (2010b), Liu et al. (2010a) and 
Shen et al. (2011) have estimated that ~0.5%-2.5% of Type 
2 AGNs at z < 0.3 are kpc-scale binary AGNs of compara- 
ble luminosities, with a relative orbital velocity > 150 km s" 1 
(see also Fu et al. 2012). Using a complementary approach 
based on a sample of 1286 AGN pairs with projected separa- 
tions < 100 kpc and velocity offsets < 600 km s" 1 , Liu et al. 
(201 1) have found a similar result of >0.5% for the frequency 
of binary AGNs on kpc scales, although the sample is biased 
against pairs with projected separations < 5 kpc. These ob- 
served frequencies based on optical surveys have been well 
reproduced by the phenomenological model presented by Yu 
et al. (201 1). The model calculates the number density of bi- 
nary AGNs based on the observed galaxy merger rate and BH- 
bulge scaling relations, under the assumption that significant 
nuclear activity is triggered only in gas-rich mergers with cen- 
tral massive black holes and only when the nuclei are roughly 
within the half-light radius of each other. 

The relatively low frequency of kpc-scale binary AGNs 
suggests that optical identification combined with X-ray con- 
firmation may be more efficient than blind X-ray searches. 
While blind hard X-ray surveys are needed to fully address 
the incompleteness of optical identification (e.g., Koss et al. 
2012), a robust determination of the frequency of kpc binary 
AGNs using existing hard X-ray surveys is still hampered by 



small sample statistics and poor angular resolution. While 
searches in the radio using eVLA and/or VLBA can probe 
much smaller scales (e.g., Tingay & Wayth 201 1; Lazio et al. 
2012), the null result from the comprehensive study of Burke- 
Spolaor (201 1) based on archival VLB I observations of 31 14 
radio-luminous AGNs may indicate that the detection yield 
could be strongly limited by the requirement of both AGNs 
being radio loud and the complication created by hot spots in 
radio jets. 

6. SUMMARY 

We have presented HST/WFC3 F336W and F105W imag- 
ing and Chandra ACIS-S 0.5-10 keV imaging spectroscopy 
of the four optically selected kpc-scale binary AGNs identi- 
fied by Liu et al. (2010a). We have further clarified the am- 
biguities concerning the nature of the ionizing sources in op- 
tically selected kpc-scale binary AGNs. We summarize our 
main findings as follows. 

1 . By combining X-ray imaging spectroscopy and star for- 
mation constraints from high-resolution U -band imag- 
ing, we have critically examined the nature of the ioniz- 
ing sources in the four optically selected kpc-scale bi- 
nary AGNs. Our new Chandra and HST observations 
confirm the binary-AGN scenario for two of the four 
targets (SDSS Jl 108+0659 and SDSS Jl 146+51 10). 
For the other two targets (SDSS Jl 131-0204 and SDSS 
J 1332+0606), the existing data are still consistent with 
the binary-AGN scenario, although the possibility of 
only one AGN ionizing both components in the mergers 
cannot be ruled out. While the new observations lend 
some further support to identifying optical kpc-scale bi- 
nary AGNs from a sample of optical candidates with 
double-peaked narrow emission lines, they also suggest 
that X-ray confirmation of optical binary candidates can 
be challenging and ambiguous, at least for heavily ab- 
sorbed sources. 

2. Combining our previous optical spectroscopy with the 
new X-ray observations, we have found tentative evi- 
dence for a systematically smaller hard-X-ray-to-[0 in] 
luminosity ratio and/or higher Compton-thick fraction 
in optically selected kpc-scale binary AGNs than in op- 
tically selected single Type 2 Seyferts. We suggest that 
the observed X-ray weak distinction may be caused by 
a combination of higher nuclear gas column (possibly 
induced by mergers) and viewing angle bias (related to 
the double-peak narrow emission line selection). 

X-ray observations of more kpc-scale binary AGN candi- 
dates are clearly needed to put our results on a firmer statis- 
tical ground. Building a much larger sample would also en- 
able the exploration of X-ray properties of kpc-scale binary 
AGNs as a function of separation and host galaxy properties. 
While deeper X-ray imaging of our targets would help better 
constrain their X-ray spectral properties, the required observa- 
tions for robust spectral modeling (i.e., > 200 counts) would 
be too expensive (likely at least ten times our exposure times 
for targets at similar redshifts with comparable [O in] lumi- 
nosities) to justify even for a small sample. Using a different 
approach, we are conducting a pilot imaging program with 
Chandra for a few kpc-scale binary AGNs at lower redshifts 
drawn from the parent sample of Liu et al. (201 1). These ob- 
servations will help better address to what extent the X-ray 
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weak tendency in optically selected kpc-scale binary AGNs is 
caused by a viewing angle bias related to double-peak narrow 
line selection (i.e., as opposed to a merger-driven gas concen- 
tration effect). 

The obscured nature of the four optically selected kpc-scale 
binary AGNs allows us to study their host galaxy stellar pop- 
ulations without much contamination from the AGN itself. In 
paper II, we will use HST WFC3 F105W and F336W imag- 
ing to characterize the detailed host galaxy morphologies and 
small-scale star formation properties. 
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